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A self-consistent spectral density approach (SDA) is ap- 
plied to the Hubbard model to investigate the possibility of 
spontaneous ferro- and antiferromagnetism. Starting point 
is a two-pole ansatz for the single-electron spectral density, 
the free parameter of which can be interpreted as energies 
and spectral weights of respective quasiparticle excitations. 
They are determined by fitting exactly calculated spectral 
moments. The resulting self-energy consists of a local and 
a non-local part. The higher correlation functions entering 
the spin-dependent local part can be expressed as functionals 
of the single-electron spectral density. Under certain condi- 
tions for the decisive model parameters (Coulomb interaction 
U, Bloch-bandwidth W , band occupation n, temperature T) 
the local part of the self-energy gives rise to a spin-dependent 
band shift, thus allowing for spontaneous band magnetism. 
As a function of temperature, second order phase transitions 
are found away from half filling, but close to half filling the 
system exhibits a tendency towards first order transitions. 
The non-local self-energy part is determined by use of proper 
two-particle spectral densities. Its main influence concerns 
a (possibly spin-dependent) narrowing of the quasiparticle 
bands with the tendency to stabilize magnetic solutions. The 
non-local self-energy part disappears in the limit of infinite di- 
mensions. We present a full evaluation of the Hubbard model 
in terms of quasiparticle densities of states, quasiparticle dis- 
persions, magnetic phase diagram, critical temperatures (Tc, 
Tn) as well as spin and particle correlation functions. Spe- 
cial attention is focused on the non-locality of the electronic 
self-energy, for which some rigorous limiting cases are worked 
out. 



I. INTRODUCTION 

Important consequences of strong correlations in nar- 
row energy bands as band magnetism, metal-insulator 
transitions or high-Tp superconductivity are normally in- 
vestigated by rather simplified theoretical models. Of 
special interest is the Hubbard modetl which describes 
itinerant electrons in a single non-degenerate energy band 
interacting via an on-site Coulomb interaction U. It in- 
corporates in the simplest way the interplay of the kinetic 
energy, the Coulomb interaction, the Pauli principle and 
the bandstructure and its consequences for the electronic 
and magnetic properties of the band system. Originally, 
the model was thought to explain the bandmagnetism of 
the transition metals and their compounds as well as the 
temperature, pressure or doping driven metal-insulator 
transition of some rtnetal oxides. The latter were first 
recognized by Motfa as contradicting conventional band 
theory (Mott transition). It is to the merit of Hubbarcfl 



to have demonstrated, how strong electron correlations 
may cause insulating behaviour of not fully occupied en- 
ergy bands. 

The in principle rather simple model provokes never- 
theless a non-trivial many body problem, that could be 
solved up to now only for some special cases .-as, e.g., 
the ground state of the one-dimensional systemQ. Other 
rigorous statements concern the possibility of collecliw; 
magnetic order in two and three-dimensional latticesBTj. 
In the first years after Hubbard's pioneering work, the 
investigation has mainly been focused on the existence 
of ferro- and_antiferromagnetic solutions of the single- 
band modclou, where, however, approximations had to 
be tolerated when tackling the respective magnetic phase 
diagram. The interest in the Hubbard model has got 
a dramatic upsurge in the recent past when its rele- 
vance tOpthe high temperature superconductivity became 
manifesfliO. It has been recognized that in the limit of 
infinite dimensions^ the electronic self-energy is wave- 
vector independenfilirtj. That simplifies but neverthe- 
less does not solve the problem—Substantial progress has 
been made by the observatior£j~L3 that for d — > oo the 
self-energy of the Hubbard model and the self-energy of 
the Anderson model are the same functionals of their re- 
spective local Green functions. However, up to now one 
does not know the exact functional form. With respect to 
the superconductivity problem, most of the recent inves- 
tigations have been done for the half filled or almost half 
filled energy band. Compared to this, relatively little ef- 
fort is concentrated on the search for magnetic solutions 
of the Hubbard model. Some interesting recent inves- 
tigations—which address the magnetic behaviour can be 
found inHhEl. 

Ferromagnetic solutions are to be expected, if at all, 
only in the strong coupling regime U/W ^irjn(W: Bloch- 
bandwidth). It has been demonstratedQ'E3 that such 
strong electron correlations are reasonable well accounted 
for by the "spectral density approach" (SDA) which basi- 
cally consists in a two-pole ansatz for the single-electron 
spectral density. The free parameters in the two-pole 
spectral density are self-consistently fixed by equating 
exactly calculated spectral moments. The advantages 
of the SDA rest on the clear and simple concept and 
its non-pcrturbative character. It turns, put to be es- 
sentially equivalent to thp_.Rpth mcthodtrEj and to the 
Mori-projector formalisnoEj. Various ^applications to 
Bose-, Fermi-, and classical systemsc3~c3 have proven 
the efficiency of the SDA. Previous applications of the 
SDA to the problem ofJhand magnetism in the Hub- 
bard model for th&,bulkE3 as well as for systems with 
reduced symmetryEj have been restricted to a local self- 
energy that could be derived self-consistently. The wave- 
vector dependent part of the self-energy, being built 



J. Magn. Magn. Mater. 170, 253 (1997) 



up by higher correlation functions referring to double 
hopping, spinflip and density correlations, has been ne- 
glected. Recentlji-jthe justification of such a neglect has 
been questionecOEJ and investigated in detail for the 
paramagnetic two-dimensional system. 

In this paper we present for the first time a full evalu- 
ation of the SDA to the Hubbard model with respect to 
the possibility of spontaneous ferro- and antiferromag- 
netisncd. Special attention is devoted to the question to 
what extent the non-locality of the electronic self-energy 
influences the stability of magnetic order. For this pur- 
pose we have organized the paper as follows. In sec. [n] we 
give a short introduction to the underlying many body 
problem. The theory of the spectral density appr oach to 
the Hubbard model is developed in detail in Sec. Ill In 
Sec.ttM we derive some analytic expressions for the case 
of half- filling (n = 1) and in the limit of strong Coulomb 
interaction U — ► oo. The results of the numerical cal- 
culations are discussed in Sec.[V| where we investigate 
para-, ferro- as well as antiferromagnetic solutions. We 
discuss the properties of the system by means of a mag- 
netic phase diagram, magnetization curves, critical tem- 
peratures, quasiparticle density of states and quasiparti- 
cle dispersions. In Sec.^3 we give a short summary and 
some concluding remarks. 



II. MANY BODY PROBLEM 

Starting point of our investigation is the single-band 
Hubbard model: 



U 



K = J2( T i] - » S ij)4* C 0<r + -j E 



H(T ' V l — (T • 



(1) 



*,J><7 



C L i c j<j) i s the creation (annihilation) operator of an elec- 
tron with spin a =|, J, in a Wannier state at lattice site 
Ri. m a — c\ a Ci a is the number operator and U denotes 
the intraatomic Coulomb matrix element. The intersite 
hopping integrals Ty are connected via Fourier transfor- 
mation to the Bloch energies: 



£ ( k ) = ^E^ e ~ 



ik^Ri-Rj 



(2) 
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Here, k is a wave-vector of the first Brillouin zone and 



Cko- the Fourier transform of Cj, 



J(+)- 



denotes the 



(anti)commutator and (...) the grand-canonical average. 
By introducing the electronic self-energy T,\ ia (E) via 
the definition 



Cka,jUY^ 



^ia^i — a 



;<L> 



e — Sk CT (-E)Gko-(-E), 



(5) 



the equation of motion of the Green function is formally 
solved: 



G^{E) 



n 



£-(e(k)- M )-£ k(7 (£) 
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The self-energy Hy L<J (E) therefore gathers all influences 
of the two particle interaction of the Hubbard model. 

The approach to the Hubbard model in this paper is 
attached to the single-electron spectral density 



SkAE) 



1 



ImG kCT (£), 
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which contains exactly the same information as the Green 
function or the self-energy. Via (inverse) photoemission 
the spectral density is directly related to the experiment. 
By use of the spectral theorem one can determine from 
the spectral density the average occupation number: 

M = {4*^)= I dEf_{E)S iia (E-ri (8) 

— OO 

+ oo 

= J dEf-(E) Pa (E). (9) 



N is the number of lattice sites. The center of gravity of 
the Bloch band is given by: 



k 



(3) 



Since the theory presented below does at this point not 
depend on the functional form of e(k) we do not spec- 
ify the underlying lattice. In sec. 5 we will discuss some 
numerical results for the bcc-lattice and the square lat- 
tice. Other lattice structures, such as sc, fee, hc-ctoo, 
fcc-<i=oo will be considered in a forthcoming papeiEj. 

All interesting single-particle properties of the system 
are determined by the retarded single-electron Green 
function: 



Here we have introduced the Fermi function f-(E) and 
the quasi-particle density of states (QDOS): 



r(E) = \s tl a{E - n) = ^- Y, S*AE - n) 



(10) 



If we assume translational symmetry, the average occu- 
pation number does not depend on the lattice site Tif. 



(riitr). 



(11) 



The total occupation number is given by n = n-j- + rij, . 

A lot of information about the spectral density can be 
drawn from its spectral decomposition (Lehmann repre- 
sentation) 
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SUE) = § 52\{E n \cl\E m )\ 2 e-^ (e<* + 1) 

x£[£-(£„-£ m )], (12) 

where 2 is the grand canonical partitition function, \E m ) 
an iV-particle, and \E n ) an (N + l)-particle eigenstate 
of the Hamiltonian. Therefore, the energy differences 
E n — E m correspond to the excitation energies required 
for adding a single electron to the TV-particle system. 

Central quantities of our procedure are the moments 
Af£™ of the spectral density Sj i<T (E), for which two equiv- 
alent representations exist. One is the usual definition via 
an energy integral 



m A = n A ]- - n A [ — n B \, - n B ^ 



-TUB- 



(17) 
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dEE n S ka (E), 
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while the other involves only commutator relations be- 
tween the construction operators and the Hamiltonian 
and is, therefore, independent of the spectral density it- 
self: 



m i n) = 

kCT N 



4£< 



ik-(Ri-Rj) 



CiaT-i 



.,H 



H,..[H,cL] 



j<?' 
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(n-p)-times p— times 

p is an integer between and n with n — 0, 1, 2, . . .. 

Besides systems with translational symmetry we also 
want to investigate antiferromagnetic structures. For 
simplicity we restrict the analytical calculations pre- 
sented in the next sections to the case of paramagnetism 
and ferromagnetism, respectively, and give only a brief 
idea of the extension to antiferromagnetism. The fol- 
lowing is meant to introduce the notation necessary to 
describe antiferromagnetic configurations. 

First we have to decompose the chemical lattice into 
two equivalent magnetic sublattices A and B. The origi- 
nal, chemical lattice can then be described by a so-called 
magnetic Bravais lattice (Rj) with a two-atomic basis 
(r Q ): 



R, 



R, 



(i = I,..., N/2; a = A,B) 



(15) 



This new labeling of the lattice sites applies, of course, 
for the creation and annihilation operators c ia<r , Ci aa as 
well. 

If we assume translational symmetry inside the mag- 
netic sublattices, expectation values do not depend on 
the lattice site R^. For example we write for the sublat- 
tice occupation number: 



(c 



(16) 



In the case of antiferromagnetic order, the dependence 
on the sublattice index a becomes important and n Aa — 
riB-a holds. This yields for the sublattice magnetization: 



Paramagnetism and ferromagnetism can still be de- 
scribed by setting n A<T = n Ba . 

All calculations are carried out in wave-vector space. 
The actual choice of the antiferromagnetic configuration, 
therefore, comes in through the single-electron energies 
e Q/3 (k) only: 



^(k) = ^E^f e " ik - (Ri " 



R, 



(18) 
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III. SPECTRAL DENSITY APPROACH 

We Mse a self-consistent spectral density approach 
(SDA)El to find an approximate solution of the Hubbard 
Hamilton operator. The method is based on a phys- 
ically motivated ansatz for the single-electron spectral 
density. The spectral density approach has proven to 
be very successful studying various many tpdy-iprob- 
lems such as Bose, Fermi and classical systemsc3~E3. The 
main advantages of this method are the physically sim- 
ple concept and the explicit non-perturbative charac- 
ter. Recent applications use the SDA successfully for 
thc_investigation of-iJic attractive Hubbard model (U < 
0)EJ, the t-J-modeO and for studying-systems with re- 
duced dimensioned, surface magnetismEj and magnetism 
of thin filmso. The SDA is believed to be a good 
starting point forr-the, investigation of high temperature 
superconductivity^ l3 . 

The work presented here, is mainly concerned with the 
magnetic properties of the Hubbard model. For the first 
time we give a complete evaluation of the SDA theory. 
The SDA can be divided into three major steps: 

(i) The crucial point is to find a physically reasonable 
ansatz for the spectral density. Some hints can be 
drawn, for example, from exactly known limiting 
cases, sum rules for the peaks of the spectral den- 
sity, spectral decompositions, etc. 

(ii) the mathematical ansatz in (i) will contain some 
free parameters, which can be fitted by the exactly 

(n) 

calculated moments M^J of the spectral density. 
All correlation functions, however, which occur in 
this procedure, have to be re-expressed by the spec- 
tral density. 

(iii) With (i) and (ii) one obtains a closed set of equa- 
tions, which can be solved self-consistently. 

If at all, ferromagnetism in the Hubbard model is to be 
expected in the strong correlation limit U > W. There- 
fore, our ansatz for the spectral density should be moti- 
vated in this limit. In the strict W = case, where no 
hopping is allowed between the lattice sites, the spectral 
density consists of two weighted ^-functions. The two 5- 
peaks are located at the energies To — /i and Tq + U — /i, 
corresponding to the excitation energies required to add 
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a (er)-electron to a lattice site where a (— <r)-electron is 
already present or not. 

The interesting question is what happens, if one allows 
for a SBiall but non-zero hopping? Following Harris and 
LangeES, four effects are to be expected: A broadening of 
the 5-peaks, a shift of the center of gravity of the peaks, 
a rearrangement of spectral weight between the peaks 
and, finally, the appearance of new, so-called " satellite" 
peaks near the energies E p .jj = Tq + p ■ U — /i (with 
p € Z \ {0, 1}). The first three effects can already be 
seen in a two site Hubbard model. Harris and Lange 
have showr£3, that the two satellite peaks at E_u and 
Ei2-u , which are closest to the main peaks, have weight 
factors of order (W/U) 4 , being therefore negligible for 
the strongly correlated system (U ^> W). The weights 
of the other satellite peaks are even smaller. Therefore, 
we can conclude that, in the case of strong correlation 
the spectral density has a two peak structure. Besides 
the Kondo-like peak— which is beyond the scope of our 
present investigationO, this decisive point in our theory 
is confirmed, for example, by recent calculations in the 
limit of infinite lattice dimensions d = odij'El 

If we neglect quasiparticle damping effects, the follow- 
ing two-pole ansatz for the single-electron spectral den- 
sity is, therefore, physically reasonable: 

S k<T {E) = hJ2 aj<r(k)*(E - E ia (k) + n). (19) 

.7=1,2 

The four free parameters Fi CT (k), F2 CT (k), ai CT (k), ck2cr(k) 
contained in the ansatz (19) can be calculated using the 
two independent representations ( |l3| ) and (JlJ) of the 
spectral moments. For this we have to calculate the first 
four moments M^J-M^J , which are given, for example, 
in ref.El. After a straightforward calculation we finally 
get: 



£(i,2)a(k) = \ 



B k - a + U + e(k) 



TV (Bk-. + U- e(k)Y + 4Un-a (e(k) - F? k _ CT ) 



aio-(k) 



E 2rT (k) - e(k) - Un_ 

£^20- — E\ G 



(20) 



(21) 



a 2CT (k) = = 1 - ai ff (k). (22) 

&2cj — ^la- 
Here, we introduced the correction term -Bk-a-, which 
contains higher correlation functions (see below). The 
correction term B k -a turns out to be of decisive impor- 
tance with respect to the possibility of spontaneous mag- 
netic order. We want to emphasize, that (E0)-(E3) repro- 
duce the exact results obtained by Harris and Langea in 
the strong correlation limit. 

It is also interesting to look at the self-energy 
^ir® A {E), which is connected to the spectral density via: 



J k CT 



S k<r {E) = h8(E-V ka {E)+n). 



(23) 



^SDA 



(E) = Un- 



E + n- B k - 



E + n-Bk-r-VQ-n-*) 



(24) 



Replacing -Bk-cr in (P4|) simply by the center of grav- 
ity T of the Bloch band, reduces Y,^ A {E) to the so- 
called Hubbard-I self-energya. It is well known that 
the Hubbard-I solution is not able to describe ferromag- 
netism. The main shortcoming of the Hubbard-I solu- 
tion with respect to the possibility of spontaneous mag- 
netic order is the fact, that the centers of gravity of the 
quasiparticle subbands (so-called Hubbard bands) in the 
QDOS are spin-independent. Allowing for such a spin- 
dependence via the correction term -Bk-o- is the decisive 
improvement of the SDA upon the Hubbard-I solution. 

The correction term -Bk-<r can be split into a k- 
independent and a k-dependent term: 



5k-f 



B_„ 



■Fk-n 



(25) 



Assuming translational invariance and hopping be- 
tween nearest neighbours only, the k-dependence can be 
separated: 

F k _ CT = (e(k) - T )F_ CT . (26) 

The remaining terms B- a and F_ CT are given by: 

1 1 ^ 

B-o-To = - -— ^T^(c l T _ (j c J _ (T (2n m - 1)), 



n_ CT (l-n_ CT ) N' 



(27) 



-[F- 



(i) 



r( 2 ) 



F 



(3)i 



, , I/ n ■ ,-, ' - n- (28) 

F_l = {rii-a-rij-a) —n^p "density correlation", (29) 

"double hopping cor.", (30) 



r — a ~ \ c ja c j-a c i-<y c i' 



F 



(3) 



\ C jo C i-{j C 3-c C i° 



"spinflip correlation". (31) 



Using (20)-(G2h we find 



The so-called spin dependent bandshift B- a may lead to 
an energy shift in the spin-spectra, allowing, therefore, 
for an exchange splitting between the spin-up and spin- 
down spectrum. 

In former investigations of the magnetic properties 
of the Hubbard model within the SDA, the bandwidth 
correction Fk_ CT was neglected and only the local part 
of the self-fej^g¥|-sKaS|-taken into account (see for ex- 
ample refs.CjSaS'LSEiluj). The bandwidth correction 
was identified to be of minor importance with respect to 
spontaneous magnetism, because the last two correlation 

(2) (3) 

functions F_£ and F_J. are effectively spin independent. 
Further, the sum of Fk-er over all wave-vectors of the 
first Brillouin zone vanishes, so no additional exchange 
splitting is to be expected. In the widely discussed case 
of infinite dimensions (d — oo) the self-energy becomes 
local, i.e., Fk-o- vanishes. Many approximations for the 
three-dimensional Hubbard model also use a local self- 
energy for simplicity. Nevertheless, in three and lower di- 
mensions the spin and wave- vector dependent bandwidth 
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correction -Fk-o- can induce substantial changes in the 
shape and the width of the QDOS. Since magnetism is 
favoured by narrow energy bands, the bandwidth correc- 
tion may also have a significant influence on the magnetic 
behaviour. One major part of this work is concerned with 
the influence of the non-locality (k-dependence) of the 
self-energy in three and lower dimensions. In particular, 
the influence of these non-local terms on the magnetic 
behaviour will be discussed. ,_. 

In one and two dimensions Beenen and EdwardsEj and 
Mchlig et. al.EZI showed for the case of paramagnetism, 
that the inclusion of the k-dependence leads to a rather 
good agreement of the quasiparticle dispersion calculated 
with the SDA and recent auantum Monte Carlo and ex- 
act diagonalization resultsOca. 

One of the more crude approximations in the ansatz 
(|19| ) for the spectral density is the fact, that quasiparti- 
cle damping is not included. The imaginary part of the 
self-energy (|24|), therefore, is identical to zero. Since the 
magnetic variables, like the magnetization to, are given 
by an energy integration over the spectral density, the 
broadening of the peaks in the spectral density should 
not be very important. However, one. has to take this 
question seriously. In a recent papeiEJ, we investigated 
the influence of quasiparticle damping on the magnetic 
stability by a proper combination of the SDA and the co- 
herent potential approximation. Guided by exact results 
in the strong coupling limit, we proposed a modified al- 
loy analogy, which allows for spontaneous ferromagnetic 
solutions. Although in the paramagnetic case the results 
obtained by the modified alloy analogy are almost identi- 
cal to the SDA-predictions, the region where the system 
exhibits ferromagnetismJ^ strongly reduced if quasiparti- 
cle damping is includedEx Further work on the influence 
of damping with respect to magnetism is in progress. 

What remains to be done in order to get a closed set 
of equations is to express the two terms B^ a and F- a by 
means of the spectral density: 

Spin-dependent bandshift: The higher correla- 
tion function (ni -c}_ .Cj_ -) which appears in the spin- 
dependent bandshift B_ a is exactly determined by the 
single-electron spectral densitycdH 



^iffW__Ci_ 



ik-(Ri-Rj) 



+oo 



k 

x JdEf_(E) [E-e(k)} S k - a (E-n). (32) 

— oo 

This leads, together with the hopping correlation func- 
tion (cj_ cr Cj_ rT ), which is given directly by the spectral 
theorem 

+ OG 

K -oo 

(33) 
to the final expression for B- a : 



U 



[E - e(k)} - 1 



S^ a (E-v)- (34) 



Bandwidth correction: The evaluation of the higher 
correlation functions contained in the bandwidth correc- 
tion F_ a turns out to be more difficult. For simplicity we 
restrict the detailed calculation to the case of the spinflip 

(3) 

correlation function F_„ . The evaluation of the density 
and the double hopping correlation functions follows ex- 
actly the same line. 



(3) 

First we rewrite F_£ as: 



F -l =-^ s iMA+±-* c i-' c i+**)i 



(35) 



where A = R; — Rj is a lattice vector which connects 
two neighbouring lattice sites. In the following we as- 
sume translational invariance and hopping only between 

(3) 

nearest neighbours. Therefore, F_J. does not depend on 
the explicit value of A. 

If we now introduce the higher spectral density A\^ 



+oo 



t'i -oo 



(36) 



Af a {t-t') = ^([ct +A _ CT c,_ (T c i+ A CT (t),cL(t')] + ) ! (37) 

it follows from the spectral theorem that F_J. is given 
by: 



F (3) = _ 



-t-oo 

^E / dEUE)A<*(E-»). 



(38) 



,(3) 



To get further information about the structure of AVI , 
we look at its spectral decomposition: 

4 3 >) - ^E< £ »iiE £ikR) 'iAV;-^ + A^) 

x(E m \ci a \E n )e~P E " [e 0E + 1) 5[E - (E n - E m )]. (39) 

The matrix elements (E n \ . . . \E m ) determine the weight 
of A^(E) at the excitation energies E = E n — E m . 

Therefore A^(E) is non-zero only if the correspond- 
ing matrix elements do not vanish. Since ( p9| ) and 
the spectral decomposition of Sk a (E) (|l^ ) involve the 
same matrix element (E m \c kcr \E n ) , the excitation ener- 
gies E n — E m are exactly the same as in S] icr (E). It is, 
therefore, reasonable to assume that the positions of the 
quasiparticle peaks in A:' (E) are also the same as in 



J. Magn. Magn. Mater. 170, 253 (1997) 



S\n a {E). However, because the other matrix element in 
( |39| ) is not equal to the one in (|12|), the distribution of 

frt\ 

spectral weight over the peaks will be different in A^ (E) 
and Sko-(E). We can, therefore, formulate an ansatz for 

4 3 J(£) = » E 0ja (k)*(£ - ^W + ^), (40) 



where Ej a (k) (j 
peaks given by ([ 



1,2) are the "old" positions of the 

and [3 CT (k) the "new" weight fac- 
tors which have to be determined. As in the case of 
Skcr(E) this can easily be done by evaluating the first 

two spectral moments of A^Ue). Note, that the same 



procedure yields the exact result (32) when applied to 
the higher correlation function {ni a c\ Cj-a) involved in 
the bandshift B- a . 

After a lengthy but straightforward calculation we find 
via (p8|) the final expression for the three higher correla- 
tion functions F ' ' . 
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F w = _ 



p(3) = _ 

1 + VQa 

For simplicity we introduced the abbreviations: 

1 



Va 



l'{\„ 



1 - n- a 

1 



(4- 



V^i — aCirrCi 



n_ CT (l 



rC,Lc 



n-a, 

r ) - n a n„ 



ja/ 



><< 



L.JCT/ 



r (l -«_£,) 



(41) 
(42) 
(43) 

(44) 
(45) 
(46) 



Since (n i ^ C! c\ IJ c : j l y 



and (c\_ a Cj-,j) are given by (32) and 
(|33|), respectively, the bandwidth correction F- a can be 
calculated by means of the sin gle- electron s p ectral den- 
sity Si«r(ip_ only. Equations @-@, @, ©-© 
and (|4l])-([46|) build a closed set of equations which can 
be solved self-consistently. 

Expressions (^l|)-(^) for the bandwidth correction are 
identicaLto the ones obtained by the method introduced 
by Rottn, both methods being, therefore, equivalent. 

For the extension of the theory to antiferromagnetic 
systems, we use an effective medium approach described 
in ref.c3. Since all further calculations are strictly along 
the lines o43 we do not give any details here. 

Paramagnetic static susceptibility: A lot of infor- 
mation about the magnetic behaviour of the system can 
be drawn from the static susceptibility: 



1 fdM\ 



(47) 



/io is the vacuum permeability, H denotes a static mag- 
netic field and 



M 



N 
V 



/i B (n T -ni) 



(48) 



is the magnetization of the system. 

For studying magnetic phase transitions we are inter- 
ested in the static susceptibility in the limit H — ► 0: 
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(49) 
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The poles of the paramagnetic static susceptibility 
X^°\n,T) indicate the instabilities of the system against 
ferromagnetic order. Therefore, all second order phase 
transitions from paramagnetism to ferromagnetism can 
be read off from \ ( n > T), which needs as an input only 
the paramagnetic solution. Besides this, the paramag- 
netic Curie temperature Q as well as the critical expo- 
nent 7 ofpthe susceptibility can easily be obtained from 
x (o)( n ,T)H 

To calculate the static susceptibility we have to include 
a Zeeman term 



H z = -HbH (rut - n i 



(50) 



in the Hamiltonian (fil) , which couples the magnetic field 
H to the local magnetic moment TOj = n^ — m^. This 
additional term does not alter or complicate the theory 
presented above. However, when differentiating the mag- 
netization with respect to the magnetic field, we have to 
take into accoimt that all expectation values are implicit 
functions of Hc5. All calculations can be done analyti- 
cally, whereas the results are rather lengthy. 



IV. LIMITING CASES 

Before we discuss the results of the numerical evalu- 
ation of the theory, we want to develop some analytical 
results for the case of exact half filling n = 1 and for 
the limit of strong Coulomb interaction U — * oo. Both 
terms, the bandshift B_ a and the bandwidth correction 
F-<j simplify considerably in these cases. Throughout 
this section we assume translational invariance, hopping 
between nearest neigbours only and a symmetric BDOS. 

First we consider the two correlation functions 
(rii-aCia c Ja) an d {4-a c i-a) which are given by (||) and 
( |33j ) respectively. The process described by (ni- a c' ilT Cj a ) 
can be understood as a hopping process of a (cr)-electron 
with the condition that one of the sites is occupied by 
a (— cr)-electron. Because of the required double occu- 
pancy this process is suppressed in the limit of large U 
and n < 1. With <ffl\ and (iSh-fEl we find indeed: 



\ c i-<j c j-ty n i<J i 



n<l,U — >oo 



0. 



(51) 



For the half-filled band, paramagnetism and T = OK, 
a simple analytical calculation involving (£32|), ([33|) and 
(P|)-(p2|) yields the relation: 
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(52) 



As U — > oo for the case of half filling the hopping cor- 
relation (cJ_g.C7_.0-) should vanish, because all hopping 
processes require an energetically unfavoured double oc- 
cupancy. In fact, (§(])- ([_§ in (§|) give for (i ^ j): 



n— 1,17— »oo 
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(53) 



We now want to focus on the spin-dependent bandshift 
B- a and the bandwidth correction _F_ CT in four different 
limiting cases: 

(i) Large U and less than half filling: With (pit) 
the bandshift _?_ CT is given by: 



B_ r 
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ZiTi 
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(cJ_ CT c,_o)+T . (54) 



__ is the number of nearest neighbours and 2_ de- 
notes the hopping integral between neighbouring lattice 
sites. It is interesting to notice that in the large U 
limit the bandshift of the (cr)-band being proportional to 
_B_- is mainly determined by the hopping of the (— cr)- 
quasiparticles. 

The density, double-hopping and spinflip correlation 
functions reduce in this limit to: 
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(ii) Paramagnetism, large U and less than half 
filling: In the case of paramagnetism we can combine 
(|5q)-(57), which yields for the bandwidth correction F-&: 
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By introducing the spin-spin correlation function (SjSj) 
this can be written as: 
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Here, the spin-spin correlation function is defined in the 
usual way: 
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(SSj) = -(S+S-) + -(Sr S +) + (S?SJ), 



(60) 



with: 
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(iii) Paramagnetism and half filling: For the exact 
half filled band (n = 1) we can use (|5-] ) to simplify the 
expressions for B~ a and F_ CT : 
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(iv) Paramagnetism, half filling and large [/: If we 

perform the large U limit starting from the exact expres- 
sion (J63J) for -F1 CT , a straightforward calculation yields: 



F_ r 
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(64) 



At this point we want to emphasize that n — ► 1 and the 
limit f — ► 00 do not commute. This can easily be seen 
with the following relation for the hopping correlation 
function which is valid for n < 1 and U — ► 00: 



(4o c j>) 
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(1-n). 



Therefore, expression ( |58| ) for F- a vanishes as n 
contrary to the correct result given by (p4|). 



(65) 
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V. DISCUSSION AND RESULTS 

The calculations in three dimensions are done for a lat- 
tice with bec structure with the tight-binding dispersion 
given by: 

e(k) = Z{F\ cos(k x a) cos(k y a) cos(k z a) + T . (66) 

a is the lattice constant, Z\ the number of nearest neigh- 
bours and T\ the hopping integral between neighbouring 
lattice sites. For the calculations in two dimensions we 
choose the dispersion of the square lattice: 



_(k) 



7 F f \ 

— - — I cos(fc x a) + cos(k y a) I + F . (67) 



Throughout this section, we fix the Bloch-bandwidth to 
W = —2Z\F\ = 2eV and the center of gravity of the 
Bloch band to To = eV. The corresponding BDOS are 
plotted in Fig.fjl. 

Out of the infinite number of possible antiferromag- 
netic structures, we consider only the relatively simple 
and most commonly used ones, AFM-(llO) and AFM- 
(AB), which are shown in Fig.pl Within the tight- 
binding approximation the hopping integrals T". read 
and the corresponding Bloch energies (|l8|) are given by: 
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A. Paramagnetism 

In Fig.pl the QDOS in three and two dimensions are 
plotted for both cases, with and without the k-dependent 
term -Fk-o- of the self-energy. The quasiparticle spectrum 
is divided into two parts, a low and a high energy region. 
These two quasiparticle subbands (so-called Hubbard- 
bands) are separated by an energy amount comparable to 
the intraatomic Coulomb matrix element U . The inclu- 
sion of the k-dependent term i^k-cr into the self-energy 
leads to a band narrowing which gets stronger with in- 
creasing band occupation. Since the centers of gravity 
of the subbands remain almost constant the energy gap 
between the subbands increases. As can be seen in Fig.||, 
the band narrowing effect gets much stronger as the di- 
mension is lowered from three to two. Since, in principle, 
magnetism is favoured by narrow energy bands one ex- 
pects the magnetic region to be enhanced by taking the 
non- locality of the self-energy into account. The band 
narrowing effect holds for the attractive Hubbardonodel 
(U < 0) as well. Contrary to the suggestions ir£3, the 
two Hubbard bands do not melt together as the non-local 
part of the self-energy is included. 

Further insight into the effect of the non-local term 
of the self-energy can be obtained from the quasipar- 
ticle dispersion which is shown in Fig.|4| for the two- 
dimensional system and moderate Coulomb interaction 
U = W. The band narrowing of the lower Hubbard- 
band is mainly due to a flattening of the dispersion near 
(1,1). The flattening extends to the saddle point (1,0) 
and halfway to (0,0). For half filling the inclusion of the 
k-dependence even leads to a local minimum of the dis- 
persion at (1,1). 

Recently, this flattening of the bands has been dis- 
cussed in j—jxmnec- 
tion with high temperature superconductivityoEil, be- 
cause similarly flat dispersions are seen in angle resolved 
photoemission experiments on hole=doped cuprate super- 
conductors. Beenen and Edwardscj show that quasipar- 
ticle dispersions like the ones in Fig.0 compare rather 
well with quantum monte carlo simulations for a 4 x 4- 
clHsteiCj'CJ. Similar results for one dimension are given 
intZI where the authors argue that the inclusion of the 
non-local terms is crucial for a good agreement of the 
model calculation with exact diagonalization results for 
a 10-site Hubbard-ring. 

For the complete information about the quasiparti- 
cle spectrum one has to look at both, the quasiparticle 
energies Ej a (k) and the corresponding spectral weights 
<x, CT (k). The spectral weight ai CT (k) of the lower quasi- 
particle subband is shown in the lower part of Fig.H. 
The spectral weight of the upper subband is given by 



CK2o-(k) = 1 — ai CT (k). If one introduces holes into the 
half filled system, there is a strong transfer of spectral 
weight from the upper to the lower subband. This weight 
transfer is particularly prominent around (1,1). 

Now, we want to focus on the bandwidth correction 
and the three correlation functions contained therein. 
In Fig. these correlation functions are plotted for the 
three-dimensional case as a function of n for three dif- 
ferent values of U. All correlation functions are negative 
and disappear, of course, in the empty band limit (n = 0) 
because of the lack of interaction partners. 

The double hopping correlation function is rather small 
and vanishes completely for U — > oo due to the un- 
favourable double occupancies. Since double hopping 
processes require a double occupied and an empty site, 
the double hopping correlation function exhibits a weak 
maximum at quarter filling and disappears at half filling 
and not too small values of U (here and in the following 
we refer to the absolute values of the correlation func- 
tions). The double hopping correlation plays, therefore, 
a minor role in the following. 

Of greater importance with respect to the k- 
dependence are the density and the spinflip correlation. 
Both correlation functions increase in a wide region of 
the band occupation n monotonously with increasing n 
and are, in this region, nearly independent of U. Close 
to half filling and for U > W, however, this behaviour 
changes and both correlation functions are strongly re- 
duced. This is because near half filling and not too small 
values of U the hopping between nearest neighbours is 
suppressed. If the lower Hubbard-band gets filled, hop- 
ping is possible only via double occupancies. 

In the case of exact half filling (n = 1), however, 
the density and spinflip correlation, and therefore -F_ CT , 
have a finite value for all U . The behaviour of n,- a (l — 
n^ a )F- a is shown in the inset in Fig.||. For U = 5.0 eV 
a smooth transition to the value at half filling is indi- 
cated, although we do not find fully converged numerical 
solutions as nS>l. For U — > oo this transition becomes 
discontinuous, indicating that the two limits n — > 1 and 
U — > oo do not commute (see Sec. IV). For not too small 
U and half filling, the numerical calculations reproduce 
F- a = —3/Zi given by (|6J) quite accurately. 

From Fig. [5] one obtains that almost all the effect of 
the density and spinflip correlation is gathered in the 
spin-spin correlation function (SjSj) which indicates an- 
tiferromagnetic correlations between nearest neighbours. 

In Fig. |6| the bandwidth correction is shown in one, two 
and three dimensions. The bandwidth correction clearly 
gets less important as the dimension increases. This fits 
perfectly well into the picture developed by the investiga- 
tions of the Hubbard model fjac infinite dimensions where 
the self-energy becomes localtll. The SDA reproduces the 
locality of the self-energy in the limit d — ► 00. 



B. Spontaneous magnetism 

Magnetic phase diagram: The global results of our 
calculations are summarized in the magnetic phase dia- 
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gram, which is shown in Fig.0. We restrict the phase 
diagram to the region < n < 1, since due to particle- 
hole symmetry the phase diagram is symmetric to the 
n = 1 axis. 

Within the SDA, we find in wide parameter regions fer- 
romagnetic and antiferromagnetic solutions besides the 
ever existing paramagnetic one. If there is more than one 
mathematical solution of the system, the physically sta- 
ble one is indicated by the minimum free energy. Sponta- 
neous magnetism is possible only for n greater than a crit- 
ical band occupation n c {U). As U increases n c decreases 
slightly without getting smaller than no = 0.54. This 
is consistent with the results obtained by Kanamorio. 
Here, we see the qualitative improvement of the SDA 
upon a simple, but .often used Hartree-Fock approxima- 
tion (Stoner model)c3, where magnetic solutions can be 
found for all band occupations, a s l ong as the Coulomb 
interaction U is strong enoughE3c3. The Hartree-Fock 
approximation is well known to overestimate the possi- 
bility of spontaneous magnetism. 

Besides a critical value for the band occupation, also 
a critical Coulomb interaction U c is needed to build up 
spontaneous ferromagnetic solutions. £/f M has a min- 
imum at n = 0.75, where, consistent with the Stoner 
criterion, the chemical potential y, lies in the region of 
maximum density of states. U c is slightly bigger 

than the Bloch-bandwidth W. However, antiferromag- 
netism in the AFM-(llO) configuration is energetically 
stable only in a small region below the ferromagnetic 
phase. Antiferromagnetism in the AFM-(AB) configura- 
tion is observed for band occupations close to half-filling. 
We find stable AFM-(AB) solutions down to very small 
values of U. However, the question, if AFM-(AB) solu- 
tions exist even for U — ► + (as for the sc lattice, see for 
exampleEJ) can not reasonably be addressed within our 
strong coupling approach. In the limit of large U , the 
region, where antiferromagnetism in the AFM-(AB) con- 
figuration is stable against ferromagnetism, is reduced to 
a small corridor around n = 1. Altogether, antiferro- 
magnetism is favoured by moderate Coulomb interaction 
U& W. 

The stability of the AFM-(AB) phase for the half filled 
band n — 1 can be explained as follows: Via second order 
perturbation theory the Hubbard model at half filling can 
be transformed into an effective Heisenberg model with 
exchange integrals equal to J^ = — 2(Ty) 2 /C7 (i ^ jjt3. 
Therefore, the total energy of the system is lowered by 
virtual hopping processes where an electron hops fraBa 
lattice site i to a neighbouring site j and back to i againt£l. 
For this virtual processes to be possible, the quasiparti- 
cles on neighbouring sites need to have reversed spins be- 
cause of the Pauli exclusion principle. Thus the lowering 
of the total energy is determined by the number of neigh- 
bouring sites occupied with reversed spin quasiparticles. 
This number is highest in the AFM-(AB) configuration, 
followed by AFM-(llO), PM and being zero for FM. This 
corresponds exactly to the results of our numerical cal- 
culations. 

Applying a modified perturbation theory (SOPT-HF: 
second order perturbation theory around the Hartree- 



Fock solution) Bulk and JelittoE3 derive a qualitatively 
very similar phase diagram. In particular, their calcu- 
lations yield similar values for the critical parameter n c 
and £7™. This is a remarkable result, because perturba- 
tion theory is usually restricted to very small values of U, 
whereas the SDA does best in the strong coupling region. 
The qualitative agreement between SDA and modified 
perturbation theory, therefore, leads to the conclusion 
that both theories produce reliable results with respect 
to magnetism for moderate Coulomb interaction U ~W 
as well. 

Ferromagnetism: In Fig. [3 the magnetization m 
of the ferromagnetic system is plotted as a function 
of n for various values of U. In addition the in- 
verse static, paramagnetic susceptibility (x (n,T)) _1 is 
shown, (x {n, T))^ 1 has either none or two roots which 
correspond to the instabilities of the system against fer- 
romagnetic order. The roots of {x^\ n , T)) -1 , which has 
been calculated from the paramagnetic solution only, in- 
dicate exactly the second order phase transitions between 
paramagnetism and ferromagnetism. As U increases the 
magnetic solutions start from lower band occupations 
and show a higher magnetization. However, there are 
two different behaviours: For U < 2.8 eV the magnetiza- 
tion curves have a maximum around n ~ 0.8 and exist 
only between the two roots of (x^H 71 : ^)) _1 - -But this 
behaviour only occurs in a parameter region, where fer- 
romagnetism is unstable against paramagnetism. Above 
U > 2.8 eV we find two ferromagnetic solutions FM1 and 
FM2, which start at the roots of (x (0) ( n > T )) -1 and ex- 
ist both until n = 1. By looking at the free energy we 
can concludeJthat FM1 is always energetically favourable 
against FM2EJ. As n increases, the FM1 solutions reach 
the ferromagnetic saturation m = n. In the limit U — ► oo 
the system is fully polarized above n — 0.65. 

In Fig. g we also show a magnetization curve for U = 
5 eV without the k-dependent term F^-a of the self- 
energy. As expected, magnetism is favoured by inclusion 
of the k-dependence. Nevertheless, in three dimensions 
the influence of the k-dependence on the magnetic prop- 
erties is rather small. The qualitative behaviour of the 
system with and without -Fk-o- is very similar. The same 
is valid in the case of antiferromagnetism. 

We now want to focus on the temperature dependence 
of the magnetization which is plotted in Fig. as a func- 
tion of n for different temperatures T and in Fig. [10|as a 
function of T and various band occupations n. As T in- 
creases the magnetization clearly decreases. Near half 
filling the two ferromagnetic solutions FM1 and FM2 
melt together. This causes first order phase transitions 
of the magnetization as a function of T for band occu- 
pations close to n = 1 (see for example the n = 0.9 
curve in Fig.JTy). Note that, again, the energetically 
stable solution is always the one with the higher mag- 
netic moment. In Fig. nn one can see that there are mag- 
netic solutions that exist for higher temperatures than all 
Curie-temperatures corresponding to second order phase 
transitions. This explains the island- like behaviour of the 
T = 600K curve in Fig.Kj. For n < 0.7 all magnetization 
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curves m(T) show a second order phase transition with a 
shape very similar to a Brillouin function, as in localized 
moment Heisenberg ferromagnets. 

From the magnetization curves m(T) as well as from 
(x^°'(n, T))" 1 we can determine the Curie temperatures 
Tc which are shown in Fig.|llj. As a function of n 
the Curie-temperatures are highest in the region around 
n ~ 0.75. Because the two ferromagnetic solutions FM1 
and FM2 melt together for very low temperatures as n 
approaches half filling, the Curie-temperatures reduce to 
zero at n = 1. In Fig.[ll](a) we also show the roots of the 
inverse static susceptibility. The system exhibits first or- 
der phase transitions as soon as both curves are disjunct. 

For a given band filling n, the Curie temperature in- 
creases with the Coulomb interaction U, but saturates in 
the limit U — > oo at a finite value which depends on n. 
We want to emphasize that this behaviour is a decisive 
improvement upon the Stoner model of band magnetism, 
for which Tc is unrealistically high and increases linearly 
with iM. 

For a better understanding of the macroscopic mag- 
netic properties like m and Tc, it is useful to look at the 
quasiparticle density of states and quasiparticle disper- 
sion. These are plotted in Fig.O for T = OK and U = 
5.0 eV. In Fig.|l|(a) the QDOS is shown for four differ- 
ent band occupations. For n = 0.55 there is no ferromag- 
netic solution and the QDOS consists of two quasiparticle 
peaks corresponding to the lower and upper Hubbard- 
bands. If the band occupation exceeds n = 0.58, how- 
ever, an additional spin-splitting occurs between p-f and 
Pl which leads to ferromagnetism. There are essentially 
two distinct correlation effects in the QDOS: An energy 
shift between the centers of gravity of the \ and the | 
subbands (Stoner shift) and a deformation of the density 
of states. The band deformation is directly connected 
to a transfer of spectral weight between the upper and 
lower quasiparticle subbands. The energy shift and the 
deformation of the band are strongly n and T dependent. 
Both effects together determine the spin-asymmetry of 
the QDOS and, therefore, the magnetic properties of the 
system. In particular, near the upper edge of the lower 
quasiparticle subband, the interplay between the band 
shift and the band deformation leads to an inverse ex- 
change splitting, i.e the J. quasiparticle states are ener- 
getically below the f states. This is because the band 
narrowing of the lower minority band overcompensates 
the Stoner shift. 

Qualitatively the band deformation and the transfer of 
spectral weight is easy to understand: The lower quasi- 
particle band results from electron hopping over lattice 
sites which are not occupied by electrons with opposite 
spin. For non-zero magnetization, the number of these 
sites becomes spin dependent. The probability of finding 
a majority electron on an otherwise empty site is en- 
hanced for m > 0. Thus the lower p^-band is narrowed 
as m increases and finally vanishes in the ferromagnetic 
saturation and n — * 1. The lower majority subband, how- 
ever, is broadened with decreasing number of minority 
particles until it becomes, for m — n, equal to the BDOS 
because of missing interaction partners. As can be seen 



analytically in the strong coupling limit, the weight of 
the lower (cr)-spin quasiparticle subband roughly scales 
with (1 — ri—a)' The weight of the upper subband, where 
the hopping of the electrons is essentially via lattice sites 
which are already occupied by a reversed spin electron, 
scales with n- a . 

The temperature dependence of the QDOS and the 
corresponding quasiparticle dispersions are shown in 
Fig.|l^(b),(c). The spin asymmetry of the partially filled 
lower quasiparticle subband decreases with increasing 
temperature. This is, again, partly due to a Stoner- 
shift and partly due to a deformation of the subbands. 
For temperatures above Tc, the spin-splitting, of course, 
vanishes, but the upper and lower parts of the spectrum 
remain separated by a gap of order U. 

The spin-splitting of the quasiparticle dispersion 
(Fig.[l2|(c)) scales with the magnetization. However, the 
actual amount of the spin-splitting strongly depends on 
the position in the Brillouin zone. There even is a region 
in k-space, where the quasiparticle energies are nearly in- 
dependent of T. Near the upper edge of the lower quasi- 
particle subband the above mentioned inverse exchange 
splitting is clearly visible. 

In the inset in Fig.|l2(c) the spectral weight of the 
lower peak in the spectral density is shown. For T = K 
the spectral weight is strongly spin dependent. As the 
temperature increases, a redistribution of spectral weight 
between the upper and lower subband leads to a decrease 
of this asymmetry. From Fig. 0(c) follows that a con- 
ventional E a (k) band-structure representation does not 
contain the full information about the underlying system, 
because details about the distribution of spectral weight 
are missing. 

The spin asymmetry of the QDOS is due to the correc- 
tion term -Bk-er in the self-energy which is given by the 
bandshift B-„ and the bandwidth correction F_ CT . With 
respect to the strong correlation limit, we have plotted 
in Fig.O the effective bandshift n^ a B^ a and the effec- 
tive bandwidth correction ri- a F- a as a function of n. 
Note that nyB^ is the effective bandshift of the lower 
t-subband. The same holds for n^ a F- a . 

For n larger than n — 0.58 the effective bandshift be- 
comes strongly spin dependent, allowing, therefore, for 
the existence of ferromagnetism. The effective bandshift 
of the majority band becomes very small and disappears 
in the ferromagnetic saturation due to n± = 0. The effec- 
tive bandshift of the minority band, however, increases 
rapidly and tends towards W/2 = 1.0 eV as n — > 1. The 
minority band gets, therefore, shifted above the Fermi 
edge. 

The spin asymmetry of the effective bandwidth correc- 
tion is less pronounced. This is because only the density 
correlation is spin dependent. Further, n- a F^ a vanishes 
in the case of half filling. This is consistent with the an- 



alytical results of section EV 



Antiferromagnetism: The sublattice magnetiza- 
tions niA for different temperatures T are plotted in 
Fig.O (a) and (b) for the AFM-(llO) and the AFM-(AB) 
configurations respectively. The most striking difference 
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to the ferromagnetic system is that the sublattice magne- 
tizations never reach the saturation vtla = n. This is be- 
cause in an antiferromagnetic system the spin up and spin 
down density of states occupy exactly the same energy re- 
gion. Filling the subbands with particles up to the Fermi 
level, therefore, always leads to a non-vanishing number 
of particles with minority spin, resulting in itia < n. 

The temperature dependence of rriA in the AFM-(llO) 
phase is quite similar to the one in the case of ferro- 
magnetism. This does not hold for the AFM-(AB) con- 
figuration, which is particularly stable near half filling. 
The two solutions AFM1 and AFM2 in the AFM-(AB) 
phase stay separated and do not melt together until fairly 
high temperatures. This can be understood by means of 
the sublattice density of states which are discussed later. 
The chemical potentials /Ui, p 2 of the AFM1 and AFM2 
solutions are separated in the AFM-(AB) phase by the 
so-called Slater gap. For AFM1 and AFM2 to melt to- 
gether, the broadening of the Fermi function has to be 
comparable to the size of the Slater gap. 

In Fig. |l5] the Neel temperatures are shown as a func- 
tion of Coulomb interaction U. In the AFM-(AB) con- 
figuration the Neel temperatures are highest for moder- 
ate Coulomb interaction, they decrease with increasing 
U in the strong coupling region and converge to a finite 
value as U —> 00. The maximum is most pronounced for 
band occupations close to half filling. As already men- 
tioned, for the half filled band the Hubbard model can 
be transformed via second order perturbation theory into 
an effective Heisenberg model with, aatiferromagnetic ex- 
change integrals J y - = — 2(Ty) 2 /[/Ello. Thus, using the 
mean field approximation for the Heisenberg model, one 
finds Xjv ~ 1/U. Close to half-filling and for not too 
small values of U the T N (U) curves in the AFM-(AB) 
phase show, at least qualitatively, this behaviour. In the 
AFM-(llO) configuration, which is built by alternating 
ferromagnetic planes, the decrease of the Neel tempera- 
tures as the Coulomb interaction increases is only weak. 
Here again, antiferromagnetism in the AFM-(llO) con- 
figuration behaves quite similar to ferromagnetism. 

The sublattice density of states of the two antiferro- 
magnetic configurations is shown in Fig.rLq. Contrary 
to the ferromagnetic system, pa^(E) and pa\{E) occupy 
exactly the same energy region due to the spin indepen- 
dence of the quasiparticle energies. The spectral weights 
which are connected to the quasiparticle energies show, 
however, a spin asymmetry leading to a spin dependent 
density of states. As for the paramagnetic and the fer- 
romagnetic system there is the splitting of the density of 
states into the lower and the higher energy region sepa- 
rated by U. In the AFM-(AB) phase an additional Slater 
gap occurs in the sublattice density of states for all pa- 
rameters, due to the high symmetry of the magnetic Bril- 
louin zone. This Slater gap is not visible in the AFM- 
(110) configuration. As the sublattice magnetization rriA 
increases, the shape of the sublattice density of states in 
the AFM-(llO) phase becomes more and more similar to 
the density of states of a two-dimensional system with 
sharp edges and a logarithmic singularity at the center. 
This quasi two-dimensionality originates from the ferro- 



magnetic ordered planes in the AFM-(llO) configuration 
(see Fig. |). 



VI. SUMMARY AND CONCLUSIONS 

We presented a self-consistent moment approach to the 
Hubbard model which is based on a two pole ansatz for 
the one electron spectral density. A complete evalua- 
tion of this theory, which does not involve any further 
approximations, leads to results identical to the method 
proposed by Roth. We believe, however, the SDA to 
be physically more transparent and better motivated, 
due to the analysis in the strong coupling limit, which 
were essential to formulate the two-pole ansatz. Since 
the SDA reproduces the band limit, the atomic limit and 
the strong coupling limit, it is expected to provide a rea- 
sonable interpolating solution covering the intermediate 
and strong coupling regime. This is confirmed in the 
paramagnetic case by the excellent agreement between 
the quasiparticle energies calculated|-with the SDA and 
recent Quantum Monte Carlo resultsEJO. For the ferro- 
magnetic case we find a phase diagram verit-similar to the 
modified perturbation theory (SOPT-HF)E^I , which does 
best for weak and intermediate Coulomb interaction. 

From the presented phase diagram we read off, that 
spontaneous ferromagnetism exists only above a critical 
value of the band occupation and the Coulomb interac- 
tion. Therefore, the SDA is much more restrictive than a 
simple mean field theory (Stoner model). Antiferromag- 
netism is most stable for intermediate coupling strength 
U ~ W and in a small region around n = 1. As a func- 
tion of temperature the magnetization curves show sec- 
ond order phase transitions away from half filling. For 
band occupations close to half filling, however, the sys- 
tem exhibits first order transitions to the paramagnetic 
phase. Note, that these first order transitions can not be 
observed by looking at the poles of the static paramag- 
netic susceptibility. The corresponding critical temper- 
atures Tc and TV lie for all parameters in a physically 
reasonable region. All magnetic properties of the Hub- 
bard model were shown to find a direct explanation in 
the (n, T, Independent QDOS. 

By investigation of the system with and without the 
non-local term in the self-energy, we find, that the in- 
clusion of the k-dependence leads to a narrowing of the 
quasiparticle subbands. The correlation functions, which 
represent the non-locality of the self-energy were dis- 
cussed in detail. Here we found, that almost all effect of 
the k-dependence is gathered in the spin-spin correlation, 
which turns out to be negative, indicating, therefore, an- 
tiferromagnetic correlation. Consistent with the d = 00 
theory, the non-local terms rapidly become less impor- 
tant as the dimension d and the number of nearest neigh- 
bours Z\ increases. For the three-dimensional bcc lattice, 
the influence of the k-dependent terms in the self-energy 
on the magnetic properties is only marginal, whereas in 
two dimensions the corrections due to the non-locality 
are essential. In two dimensions, however, we did restrict 
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our investigations to the paramagnetic case, since spon- 
taneous magnetic order in the two-dimensional Hubbaxd 
model is excluded by the Mermin Wagner theoremE-m. 
Very recent calculations indicate, that for the three- 
dimensional sc lattice the influence of the non-local terms 
with respect to ferromagnetism is much more important 
than for the bcc lattice. We will discuss this fact in a 
forthcoming paper, where we investigate the influence of 
the lattice structure on ferromagnetism in the Hubbard 
modeO. 
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FIG. 1. Bloch density of states (BDOS) for (a) a three-dimensional bcc lattice and (b) a two-dimensional square lattice. 
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FIG. 2. Decomposition of the bcc lattice into two equivalent sublattices. The AFM-(llO) configuration is built in such a 
way, that the (HO)-planes alternately belong to the sublattice A and B. In the AFM-(AB) configuration all nearest neighbours 
of a given lattice site belong to the other sublattice. The vectors ai,2,3 are the translation vectors of the respective magnetic 
Bravais lattice and a is the lattice constant. 
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FIG. 3. Paramagnetic quasiparticle density of states (QDOS) as a function of the energy E in (a) three and (b) two dimensions 
for different band occupations n and the Coulomb interaction U — 5.0 eV. Solutions with (solid lines) and without (broken lines) 
the non-local part of the electronic self-energy are shown. The vertical lines indicate the positions of the chemical potential jj,. 
(W = 2.0eV, T = 0K) 
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FIG. 6. Bandwidth correction F_ CT in three, two and one dimension as a function of the band occupation n. The values of 
F-„ for half filling n = 1 are given on the right. (W = 2.0 eV, T = OK) 
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FIG. 7. Magnetic phase diagram for the three-dimensional Hubbard-model with bcc lattice structure, n is the band occu- 
pation and U the intraatomic Coulomb matrix element. The bandwidth of the BDOS is fixed to W — 2.0 eV. Four phases 
are considered: Paramagnetism (PM), ferromagnetism (FM) and antiferromagnetism in two configurations AFM-(llO) and 
AFM-(AB). (T = 0K) 




FIG. 8. Magnetization m as a function of the band occupation n for different Coulomb interactions U. In addition, the 
inverse static paramagnetic susceptibility (x '(n,T))~ is shown in arbitrary units. (W — 2.0 eV, T — OK) 
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